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Abstract. This paper concerns the restricted 3-body problem. By applying topological meth- 
ods we give a computer assisted proof of the existence of some classes of periodic orbits, the 
existence of symbolic dynamics and we give a rigorous lower estimate for the topological entropy. 

1. Introduction 

The problem of both restricted and full A'^-body systems has such a long story that it is 
impossible to give an extensive bibliography here; we refer the reader to the classical texts 
[MH, M, SM, S]. 

This paper concerns the study of periodic and chaotic solutions of the planar restricted 3-body 
problem. If we assume that the primaries orbit around each other with period 27r and we use a 
rotating reference frame, i.e. if we use synodical coordinates, then the motion of the third body 
is described by the following system of second order differential equations: 

X + 2y = 
y 2^ — 5 



(1.1) 

where 



n{x, y) = — + !^+ + + a 



nil and m2 are the masses of the primaries, C is an arbitrary constant and R2 depend on 
mi and m2- 

The methods we employ are topological with computer assistance. This kind of computer 
assisted proof has been introduced in [MMl, MM2, MM3], while the topological methods em- 
ployed here have been introduced in [Zl, Z2], see also the applications in [AZ, GZ, Z3]. We point 
out that, although some of the proofs we give require computer assistance, they are rigorous 
and they can be easily reproduced on any recent computer. To this purpose, in the last section 
we provide all necessary information, mostly by referring to other papers, and we also provide 
on the web a Mathematica version of the algorithms employed. 

The main result presented in this paper is the proof of symbolic dynamics on fourteen symbols 
which corresponds to the existence of orbits, periodic or non-periodic, that come close to five 
different periodic orbits in any prescribed order. This result yields a lower estimate for the 
topological entropy of the system, and it is therefore a proof of its chaotic behavior. To the 
author knowledge, no rigorous estimate of the topological entropy for the Poincare map of 
this system is available in the literature. This paper serves other purposes as well. We show 
how to extend the computer assisted techniques developed in the papers cited above to the 
planar restricted 3-body problem, which presents different features and difficulties and requires 
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new techniques; to this purpose we introduce some new topological tools and computational 
methods. We give a rigorous proof of the existence of a class of periodic orbits at different 
energy levels: the existence of such periodic orbits (and many more) is well known, but to the 
author knowledge the proof is sometimes purely numerical, with no mathematical rigour, or 
perturbative, with no estimate on the range of validity of the perturbation parameter. We also 
provide a very narrow estimate of the locations of the intersection of such periodic orbits with 
the line connecting the primaries. In summary, this paper provides both the methods and some 
examples on how the computer assisted techniques can yield results for the planar restricted 3- 
body problem. We point out that these results are not of pertubative nature, and to the author 
knowledge the result on symbolic dynamic is not accessible by purely analytical methods; in 
particular the problem is not treatable by Melnikov's method. 

A result on chaos for the planar restricted 3-body problem has been recently presented in 
[SK] with a different method, also based on computer assisted techniques, but the authors do 
not claim to give a rigorous proof. Indeed, they use, in their words, "realistic estimated upper 
bounds" for the errors made in the numerical computation of the Poincare map. Such method 
found a rigorous application on a different system in [KMS]. Here we estimate rigorously all 
computational errors and we do provide rigorous proofs of all the theorems we state. 

The paper is organized as follows: in Section 2 we provide a brief introduction to the problem; 
in Section 3 we introduce the Poincare maps that are used in the proofs and explain a symmetry 
of the system which is widely employed in the proofs; in Section 4 we explain our method 
for detecting and proving periodic orbits and present the results we obtain. Section 5 is the 
main portion of the paper. Here we introduce the topological and computational methods and 
we describe the results on symbolic dynamics and on topological entropy. Ideas on further 
developments are given in Section 6 and details on the computer assisted proofs are in Section 
7. 



2. Description of the system 

We first derive briefly equation (1.1). It is well known that the two body problem admits a 
solution where both bodies move in a circular counterclockwise motion around their center of 
mass (and of course it also admits the symmetric clockwise solution). We call the two bodies 
Pi and P2 (the primaries). If Pi and P2 have mass mi and 7712, then the circular solution with 
minimal period 27r minimizes the Lagrangian functional 

r-27r 



f( \ ■ |2 , "^2| . |2 , 

Jo 2 2 \xi 



mim2 

—at 



X2\ 

defined on i?^g^([0, 27r]), the space of 27r-periodic functions whose weak derivative is square- 
integrable. 
Take 

xi = Pi (cost, — sint) and X2 = R2{ — cost, sint); 
in order to compute Pi and P2 we have to minimize 

(2.1) = + + 

Setting as the center of mass we have P2 = ^Ri, and by minimizing /(Pi,P2) we get 
Ri = .^."1^2/3 and P2 



(mi +9712)^/^ (mi+m2)^/^' 
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The Lagrangian for the restricted 3-body problem is L{w) = — V{w) where 

V{w) - 



\W — Xi\ \W — X2\ 

It is convenient to use synodical coordinates, i.e. a reference frame where the primaries sit still. 

Set w = Rv, where R = ( "^^f * sint \ let J = | ? ^"'^ | be the standard symplectic 
' \^ -smi cosi y \1 J ^ ^ 

matrix; then w = Rv + Rv = R{v — Jv), therefore the Lagrangian is given by 

L{v) = ^\v\'' -iv,Jv) + n{v), 

where 

It is well known that H{x, x, y, y) = aP' + — 2^{x^ y) is an integral of the motion (the Jacobi 
integral), therefore the motion takes place on the manifold H{x,x,y,y) = h. With some abuse 
of language, we call the Jacobi integral H{v) = |up — 2Q{v) the energy. The Euler-Lagrange 
equations are: 

(2.3) i-2Ji)- VVl{v) = 0. 

This system admits five equilibrium points. Three such points are aligned with the primaries 
and are called Li, L2 and (the collinear equilibrium solutions). One of the collinear points 
lies between the primaries; let that point be Li. All the collinear points are saddles for Jl. The 
remaining two equilibrium points are called L4 and L5; they are the absolute minima of O and 
they both sit at the vertices of an equilateral triangle with the primaries at the other vertices 
(the equilateral equilibrium solutions). 

We only consider the case mi = m2- The reason for this choice is that this case has been 
already widely considered and many families of periodic orbits for this case are known to ex- 
ist. Nonetheless, rigorous estimates on the location of the periodic orbits, results on symbolic 
dynamics and topological entropy and even a rigorous proof of existence of the periodic orbits 
are new to the author knowledge. Furthermore the equal masses case is, in some sense, the 
opposite of a "perturbative" case with one primary much less massive than the other one. In 
this case Ri = R2 = 2^^/'^ and we choose the constant C = —2^1^ in order to have ^{Li) = 0. 
By a direct computation it is easy to see that there exists Hq = —2i}{L2) = — 2ri(L3) ~ .8623 
such that the region {v : 2Q{v) + h > 0} where the motion can take place is split in exactly 
two connected regions, one bounded and the other unbounded, if and only ii < h < Hq. We 
only consider values of h in the range [0, ho) and we look for trajectories in the bounded region. 
Note that, if /i = 0, the bounded region is equal to the union of two closed sets whose only 
intersection is the origin and no trajectory can intersect both sets; if < /i < /iq the bounded 
region is homeomorphic to a ball (if we include the singularities in the region); if /i < there are 
three connected regions and finally if h > ho the region {v : 2^l{v) + h > 0} is connected, but 
not simply connected. The following pictures represent the curves {2^l{v) + h = 0} for different 
values of h (the small disks represent the primaries): 
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h = .3 h = .9 



3. The Poincare maps 

In order to study the system at some fixed energy h we consider the Poincare return map 

P : D{P) CM^ defined in the fonowing way. Given {x,Px) such tliat x ^ {Ri, ~ R2} and 

2il(x, 0) + h-pl > there exists a unique positive value of Py — Py{x , Px ) = ^/2Q{x,0) + h - pi 
such that H{x,px,0,Py) = h. Let (p{x,px; ) : M — > M'^ be the solution of the equation (1.1) 
with initial conditions 2;(0) = x, x{0) = px, y{0) = 0, y(0) = py and call ipi, i = 1, ... ,4 its 
components. By definition ip:^{x.px;0) = and ip4(x,px;0) > 0, therefore (p3{x,px;t) > for all 
positive and small t. If there exists a time Ti such that (p3{x,px;Ti) = and ip3{x,px;t) > 
for all t G (0,Ti), and a time T2 > Ti such that ips{x,px;T2) = and (p3{x,px;t) < for all 
t G (Ti,T2), then we define 

P{x,Px) = {m{x,Px;T2),ip2{x,Px;T2)). 
In other words, P is the return map on the section y = 0, py > 0. 
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In fact, for different reasons which we point out later, we find it useful to consider the half 
Poincare maps Hi : D{Hi) CR^ and H2 : D{H2) CR^ ^R^. Hi is defined by 

Hi{x,px) = {ipi{x,px;Ti),ip2{x,px;Ti)), 

where the maps ip and the time Ti are as before, while H2 is defined by 

H2{x,p^) = {ipi{x,p^;Ti),(p2{x,p^;Ti)), 

where <p{x,px; •) : M — > is the solution of the equation (1.1) with initial conditions a;(0) = x, 
±{0) = px, y(0) = 0, y(0) = -^y2n{x,0) + h-pl and Ti > is such that ip^{x,pr,-Ti) = 
and (^3(2;, Pa;; i) < for all t G (0,Ti). Of course, D{P) C D{Hi), Hi{D{P)) C D{H2) and 
P = H2 o Hi. In the following we refer to the maps Hi and H2 as the first and second half 
Poincare maps respectively. In other words the map Hi (resp. H2) is the transition map from 
the section y = 0, py > to the section y = 0, py < (resp. from the section y = 0, py < to 
the section y = 0, py > 0). 

By inspection it is easy to see that, if {x{t),y{t)) is a solution of (1.1), then {x{t),y{t)) : = 
{x{—t),—y{—t)) is also a solution. This implies that Hi{xi,pi) = {x2,p2) is equivalent to 
H2{x2, -P2) = {xi, -pi) and P{xi,pi) = (2:2,^2) is equivalent to P{x2, -P2) = {xi, -pi). 



4. Periodic orbits 

A standard method for studying periodic orbits consists in looking for fixed points of the 
Poincare map P. By the symmetry of the system considered at the end of the previous section 
we infer that i7i(a;i,0) = (3^2, 0) yields i72(2;2,0) = (xi,0), which in turn implies that (a;i,0) is 
a fixed point for the Poincare map. On the other hand, if the system admits a periodic orbit 
which crosses orthogonally the x-axis at some point xi, then, by definition of Poincare map, 
P(a;i, 0) = {xi, 0), and this is possible only ii Hi{xi, 0) = {x2, 0) for some X2. It turns out that it 
is also useful to consider whether there exist points xi and X2 such that P{xi, 0) = (x2, 0), with 
xi ^ X2. By the same reason as before, these points correspond to periodic points of period 2 for 
the Poincare map, hence to periodic trajectories for the system crossing the y = hyperplane 
at two different points in each direction, see the pictures of the orbits Di and D2 below. 

Let f{x) be the second component of i?i(a;,0) and g{x) be the second component of P(a;,0). 
In order to find fixed or period 2 points for the Poincare map we can look for zeros of the 
function / or g. We remark that, even by considering the map g only, one can still find all the 
fixed points for the Poincare map. On the other hand it may happen that P{xi,Q) = (0:2, 0) 
with xi very close to X2, in the sense that \xi — X2\ is smaller than the numerical error. In 
this case it is impossible to find out whether xi is a fixed point or a periodic point of period 
2, without considering either the derivative of the Poincare map or the half Poincare map. For 
this reason it is convenient to study both the map / and the map g: first we look for zeros of 
/, i.e. fixed points of P, then we can look for zeros of g which are not zeros for /, such points 
yielding periodic points of minimal period 2 of P. The following picture displays a numerical 
computation of f{x) with x G ( — .62, .62) and h = .1: 
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The picture strongly suggests that the Poincare map has three fixed points in the set {{x, 0) : 
X € (— .62,-62)}, which is exactly the result we show rigorously in the next section, where we 
also give a narrow bound on the position of such points. Our strategy to detect and prove 
the existence of periodic orbits is as follows: first we choose some values for h and compute an 
approximate image of the x axis through the maps / and g as in the previous picture; in this way 
we can spot the places where the intersections should be located. Assume that the numerically 
computed graph of / crosses the axis in a neighborhood of some point x: we conjecture the 
existence of a fixed point for P nearby some point {x, 0). To prove the conjecture we choose xi 
and X2 such that xi < x < X2 and both xi and X2 are very close to x. Then we compute the 
rigorous half Poincare map Hi at (xi, 0) and (x2, 0). If we can prove that the second component 
of Hi{xi,0) has opposite sign with respect to the second component Hi{x2,0) and that the 
segment joining the two points belongs entirely to the domain of Hi, then by the continuity 
of the half Poincare map we have proved that there exists at least a point xi < x < X2 such 
that Hi{x,0) lies on the x axis, therefore a periodic orbit intersects the a:-axis orthogonally at 
{x, 0). On the other hand, if we prove that some portion of the a;-axis is mapped away from the 
X-axis itself, then we have a proof that there are no periodic solutions which cross the x-axis 
orthogonally in that section. Then we search for points of period 2 for the map P, that is we 
study the numerically computed graph of g and look for intersections with the a;— axis. If we 
find two points xi, X2 such that the second components of P{xi,0) and P{x2,0) lie on the 
opposite sides of the x-axis, the set [a;i,a;2] x {0} belongs to the domain of P and P{[xi,X2],0) 
does not intersect the set [a;i,a;2] x {0}, then we have a proof that there exists at least a point 
xi < X < X2 such that x is the intersection of an orbit of minimal period 2 with the a;-axis. 

The case we consider is usually referred to as "the Copenhagen orbits", from the results of the 
Observatory of Copenhagen (see [S] and references therein). We recall that, by our choice of the 
constant C in (2.2), /i = is the energy of the stationary solution at the Lagrangian point Li (the 
origin). Since the problem has the additional symmetry consisting in switching the primaries, 
we only look for orbits that cross the portion of the x-axis between the primaries with positive 
speed in the y direction. The lowest value of h we consider is 0, when the admissible region is 
split in two parts touching at the origin and no trajectory can enter both regions. The highest 
value of /i we consider is 0.8, since at slightly larger value the bounded admissible region touches 
the unbounded part and some trajectory starting close to the primaries may be unbounded. 
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h 

.1 
.2 
.22 
.24 
.26 
.28 
.3 
.4 
.5 
.6 
.8 



Si 
0.3158 
0.3028 
0.2883 
0.2851 
0.2818 
0.2784 
0.2749 
0.2712 
0.2443 



0.4399 

0.4365 
0.4330 
0.4323 
0.4316 
0.4309 
0.4301 
0.4294 
0.4257 
0.4218 
0.4178 
0.4093 



S2 



0.02697 
0.03894 
0.04102 
0.04303 
0.04497 
0.04687 
0.04873 
0.05751 
0.06575 
0.07369 
0.08922 



L 



0.04423 
0.04600 
0.04775 
0.05615 
0.06412 
0.07185 
0.08716 



0.04456 
0.04651 
0.0484 
0.05730 
0.06550 
0.07345 



D2 



0.06963 
0.06737 
0.06689 
0.06712 

0.07199 
0.07870 
0.08595 



0.1043 
0.1205 
0.1343 
0.1470 
0.2085 



U2 



Table 1 



The main result of this section is that the system (1.1) admits periodic solutions as shown in 
Table 1. In the left hand side column the energy level is displayed, while the remaining columns 
represent the a;-coordinate ±5 • 10^^ of the intersection of the orbits with the portion of the 
a;-axis between the primaries (the L, Di and D2 orbits have two such intersections; for the L 
orbits we only consider the intersection with positive y velocity, the other being symmetric with 
respect to the origin; for the Di and D2 orbits we consider the only orthogonal intersection). 
More precisely, the following theorem holds: 

Theorem 4.1. For all energy values displayed in column h of Table 1, the system (1.1) admits 
at least a periodic solution for each value printed in the remaining columns. Such solution 

crosses the x-axis orthogonally twice. The x-coordinate of the intersection with positive y-velocity 
lies in the interval centered in the position given in the table with width 10~^. The solutions 
corresponding to values in the columns Si, S2, L, Ui and U2 do not intersect the x-axis at any 
other point, while the solutions corresponding to values in the columns Di and D2 intersect the 
X-axis twice at another point with the same negative x-velocity and with opposite y-velocity. The 
orbits in the column S2 are retrograde around Pi; the orbits in the classes Si, Ui and U2 are 
direct around P2. The orbits in the class L are retrograde around Li. 

Proof. Orbits Si, S2, L, Ui and U2. We checked by computer assistance (see Section 7) that 
the second components oi Hi{x + 5 ■ 10~^, 0) and Hi{x — 5 ■ 10~^, 0) have different sign and the 
interval set [x - 5 ■ W-'^,x + 5 • lO""^] x {0} is in the domain of Hi, where x is any of the values 
in the columns Si, S2, L, Ui and U2. According to the argument presented at the beginning of 
this section, this suffices to prove the existence of a fixed point of the Poincare map P and hence 
the existence of a periodic orbit. To prove that an orbit in the column 52 is retrograde around 
Pi we proceed as follows. We compute the trajectory of the set [x — 5 ■ 10^^, a: + 5 • 10^^] x {0} 
until it crosses the Poincare section and we check that the angular velocity with respect to Pi is 
always strictly positive. This also implies that the a;-axis is crossed only once in each direction, 
and by construction and the symmetry of the system the intersections are othogonal. To prove 
that the orbits in the classes Si, Ui and U2 are direct around P2, the orbits in the class L are 
retrograde around Li and they also cross the a;-axis only once in each direction we follow an 
analogous procedure. 
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Orbits Di and D2. Same argument with the map P. To check that the periodic points 
of the Poincare map corresponding to these orbits have minimal period 2 we checked that 
P{[x - 5 ■ 10-^ a; + 5 • 10-''],0) is mapped away from [a; - 5 • 10-^a; + 5 • IQ-''] x {0}. The 
symmetry of the system forces the two non-orthogonal intersections with the a;-axis to occur at 
the same point, with the same rr- velocity and with opposite y- velocity. □ 

Remark 4.1. Orbits belonging to the same column appear to belong to the same class, in the 
sense that they have the same linear stability and they have the same winding number with respect 
to the primaries and the Lagrangian point Li . More precisely, all orbits are unstable, except for 
orbits Si and S2 which are linearly stable. The result on the stability is purely numerical. 

The orbits in the class L are well known, their trajectory is very close to an ellipse with center 
at the Lagrangian point Li. They branch out from Li, in the sense that as /i — > they collapse 
to Li. These are the only orbits considered in this paper which do not wind around any primary, 
while they wind around Li. The orbits in the class Di and D2 correspond to points of period 
2 of the Poincare map. In the following pictures the orbits Di, D2 and L at energy level /i = .5 
are displayed (the small disk represents Pi). 

y y y 




The same orbits are plotted together in Figure 1: note that part of the trajectories of both 
Di and D2 are very close to the trajectory of L, giving a visual image of the strong instability 
of the system. 

Remark 4.2. The existence of unstable orbits around the other Lagrangian points is also well- 
known, but they occur only for higher values of the energy; we do not investigate those orbits, 
but we claim that the methods presented in this paper could be used to study those orbits and, 
possibly, chaotic dynamics involving those orbits as well. 

In Figure 2 we plot all periodic orbits at energy level h = .3: note that the shape of the orbits 
Di, D2 and L are quite similar to the same orbits in Figure 1, supporting the claim that they 
belong to the same branch. 

5. Chaotic dynamics 

We apply the method developed in [Zl, Z2, Z3], see also [AZ] where the Henon-Heiles Hamil- 
tonian was concerned. We only consider the system at energy h = .3. We prove the existence 
of symbolic dynamics on 14 symbols and we give a lower estimate for the topological entropy. 




The computation of rigorous estimates in the (restricted) 3-body problem presents more 
difficulties than the cases previously considered. The first is due to the fact that the Poincare 
map is not defined in a connected region, indeed we have to exclude at least the lines x = Ri 
and X = — i?2- This could be avoided by using some kind of regularization, but we prefer to 
keep the original coordinates, both because they are far more intuitive and because the Levi- 
Civita transform (or its variations) is not one to one. Furthermore, even after some kind of 
regularization, it is far from trivial to determine what is the domain of the Poincare map. But 
the main problem is due to the fact that, even using a sofisticated algorithm to compute the 
rigorous bounds, such bounds turn out to be very large, particularly in the Px direction. The 
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reason for such large bounds are not only the inevitable computational errors, which in theory 
can be as small as we like (by taking a smaller time step at the price of a slower computation), 
but particularly the wrapping effect. See [AZ, GZ] and the references given there for a discussion 
on this topic. Here we just want to point out that although it is possible by a change of variable 
to get rid of the singularities, it is not possible to get rid of the wrapping effect in the same way, 
at least not with the same change of variable. 

The first trick we had to adopt consists in computing only half of the Poincare map at a 
time, and this is the second reason, in fact the most important, for introducing the maps Hi 
and H2. Indeed the wrapping effect is usually exponential, therefore by considering about half 
trajectory time it is rather drastically reduced. By this method we obtain a great reduction 
of the error: of course we have to pay the price of a larger number of computations, but the 
trade-off is very positive. The second trick that proved to be essential consists in computing the 
inverse of the map instead of the actual map for some checks, see the definition of back-covering 
below. This is important whenever some trajectory starts away from the primaries but comes 
close to one of them at the intersection with the Poincare plane, indeed in such cases a very 
short time step is necessary to keep the error small, since close to the singularities the speed 
undergoes a strong variation while the particle crosses the Poincare plane. The particle is in 
fact a point, but in order to compute the Poincare map we have to consider the envelope of its 
position with the error bounds, hence it takes a finite time to cross the plane. It turns out that 
the inverse trajectory, starting close to a primary and ending away from both primaries, raises 
a much lower error. 

5.1. Topological tools. Definitions 5.1 and 5.3 were introduced in [AZ]. Definitions 5.4 and 
5.5 are introduced here for the first time to deal with the new difficulties of the restricted 3-body 
problem. Lemma 5.3 and Corollaries 5.4 and 5.5 are also new. 

Definition 5.1. A triple set (or t-set) is a triple N = (|iV|, iV^ iV) of dosed subsets of 
satisfying the following properties: 

1: I TV I is a closed parallelogram 

2: and iV" are closed half-planes 

3: the sets N^^ := fl \N\ and N^*^ := iV fl \N\ are two nonadjacent edges of \N\ 
We call \N\, N\ , N^'^' and N^'^ the support, the left side, the right side, the left edge and 
the right edge of the t-set N respectively. One observes that \ (|7V| Li U-ZV) consists of two 
disjoint sets. We call and N'' the closure of such sets (the top and bottom sides of the triple 
set). We also set N^" := iV* \N\ and N''"' := fl \N\ (the top and bottom edges of N). 

Remark 5.1. Since all theorems we use concern topological properties of the t-sets, then one 
can choose any continuous deformation of the t-set defined above, obtaining the same results. 
An example of the sets we actually use, whose exact definition is given in Section 7, is given in 
Figure 3. 

To exploit the symmetry of the system we give the following definition: 

Definition 5.2. Let M be a t-set. We define its symmetric image with respect to the x-axis 
M as follows: if S : -5- is the map defined by S{x,y) = {x,-y), let \M\ = S{\M\), 
M' = S{M^), M*" = S{M^), and the remaining definitions follow as in Definition 5.1. 

Definition 5.3. Let / : f2 C ^ be a map and let Ni and N2 be two triple sets. We say 

that Ni f- covers N2 (Ni =^ N2) if: 
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Figure 3. An example of a t-set 



a: /(|iVi|) C mt(A^^ U \N2\UN^) 
b: either f{N[^) C int(iV^) and /(A^f) C mt{N^) 
or f{N['=) C mt{N^) and /(7V[^) C int(7V^) 

The following lemma says that we can reduce the condition (a) in the above definition to the 
boundary of |7Vi| if we know that the map / is defined on |A"i| and it is injective. This lemma 
plays a very important role in the computer assisted verification of the covering relations, as it 
reduces the computations to the boundary of |A^i| (see Section 6 in [GZ], for more details). 

Lemma 5.1. Let / : C — )■ be a map and let Ni and N2 be two triple sets. Assume that 

f is an injective map on \Ni\, then Ni =^ N2 if and only if 

a': f{d\Ni\) C int(A'^ U jiVal U N^) 
b: either f {N[^) C int(iV^) and f{Nl^) C int(A'2^) 
or f{N[^) C int(JV2^) and f {Nl^) C int(7V^) 

As we pointed out above, in some circumstances it is easier to compute the inverse flow, then 
the direct flow. Furthermore, we need to deal with the half Poincare maps. To this purpose we 
give the following definitions: 

Definition 5.4. Let Ni and N2 be two triple sets. We say that Ni f — backcovers N2 (Ni 

■^= N2) whenever: 

a: / : C — >■ C is a homeomorphism 
b: \N2\ C n2 

c: /-^(a|A^2|) C mt{Nl U |A^i| U N^) 
d: either f'^iN^") C int(Arf) and f-^{NP/) C mt{N^) 
or f-^iNl") C int(A^I') and /"^(A^l^) C int(A^|) 

Definition 5.5. Let Ni and N2 be two triple sets. We say that Ni generically f — covers N2 
(N\ <;=^ N2) if Ni f — covers N2 or there exists n > I t-sets Mi, « = 1, . . . , n and n + 1 maps 
Qj, j = 0, . . . ,n such that f = Qj o ■ ■ ■ o Qi o go, Ni =^ Mi, Mi =^ Mj+i for all i = 1, . . . , n — 1 
andMn ^N2. 



12 



G. ARIOLI 



Remark 5.2. Although covering, backcovering and generic covering occur often simultaneously 
and they are indeed a very similar phenomenon, they are not equivalent. In fact it may even 
happen that a map f is not defined on the whole support of Ni, and still Ni f—backcovers N2 
or Ni generically f— covers N2 ■ On the other hand, the result of a backcovering or of a generic 
covering relation is very similar to the result of a covering relation as far as the results we are 
interested in are concerned, see Corollaries 5.4 and 5.5. 

The following theorem follows immediately from Theorem 4 in [Z3] and together with Theorem 
5.6 it is the main topological tool we use: 

Theorem 5.2. Given n t-sets Mi C and n continuous maps fi'.Mi^'^, such that 

Mq Ml M2 ^ M2 . . Mq = Mr,, 

then there exists x G int|Mo|, such that f/^o ■ ■ ■ o fio fo{x) G int|Mjt+i|, for A; = 0, . . . , n — 1 and 
X = fn-l o ■ ■ ■ o /i o fo{x). 

By the following lemma we can extend the previous theorem to the generic backcovering, see 
Corollary 5.5. 

Lemma 5.3. Let Oi and O.2 be two open sets o/M^, let / : fii — )■ be a homeomorphism and 
let Ni and N2 be two triple sets such that \Ni\ C Oj, i = 1,2. If Ni f—backcovers N2 then there 
exists a t-set K C such that Ni id— covers K, and K f— covers N2 (id is the identity map 
in 

Proof. By definition of backcovering, if TVi /— backcovers N2, then /^^(|A^2|) looks as the light 
grey rectangle in the Figure 4, therefore it is possible to define a t-set K such that the boundary 
of its support is very close to the boundary of /^^(|A^2|)5 in such a way that 7Vi2c?— covers K and 
K /—covers N2. To prove this, just consider a set K' as the dark grey rectangle in the picture 
or, more precisely, if a and h are the length of the sides of |iV2|, let K' be the parallelogram with 
the same center as |iV2|, the sides parallel to the sides of |A^2| and such that the length of its 
sides is (1 — e)a and (1 + e)6, < £ < 1. If e is small enough, K' C ^^2 and f^^{K') is defined. 
Call the edges of K' top, bottom, left and right according to the corresponding edge of N2. 
Now define K as follows. The edges are the counterimage of the edges of K' . By construction, 
definition of backcovering, compactness of K' and continuity of /~^, if e is small enough the 
top (resp. bottom) edge of K lies in the top (resp. bottom) side of A^i and both the right and 
the left edges are in int(A''| U \Ni\ U iVf), hence it is possible to define the left (resp. right) side 
of K in such a way that A^^i id— covers K. By construction K /—covers N2 and the proof is 
complete. □ 

Corollary 5.4. Let Ni and N2 be two triple sets and let f be a homeomorphism. If Ni 
f — backcovers N2, then Ni generically f — covers N2. 

Corollary 5.5. Given t-sets Mi C and n continuous maps fi : ilu C ^ ^21 C such 
that 

Mo <^ Ml 4^ M2 <^ M2 . . . ^ Mo = Mn, 
then there exists x G int|Mo|, such that fk°'"°fi°fo{x) G int|Mfc_|_i|, for k = 0, . . . ,n — 1 and 
X = fn-i o • • • o /i o /o(a;). 

Assume that we have n t-sets A^j, « = 1, . . . , ra, with some covering relations. Let N = \J-\Ni\; 
the following definitions are standard. 
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Figure 4. In light grey / ^(iV2), in dark grey K' = f[K) 

Definition 5.6. Let f be injective. The invariant set of N is defined by lBY{N,f) := {x & N : 
f{x) G TV for all i G Z}. 

Definition 5.7. The transition matrix T{j,i), i, j = 1, . . . ,n, is defined as follows: 

I otherwise 

Let E„ be the set of bi-infinite sequences of n symbols 

Definition 5.8. A sequence {x^} G S„ is said to be admissible ifT{xk+i,Xk) = 1 for all k. We 
denote by C S„ the set of all admissible sequences. 

Definition 5.9. Assume \Ni\ fl |A''j| = 0, for i / j. The projection tt : lnv{N,f) — )■ Ha is 
defined by setting Tr{x)i = j where j satisfies f^{x) G |7Vj| for all 2 G Z. 

The set inherits the topology from E„; the shift map a : T,a is continuous. We 

prove a semiconjugacy between a and /, i.e. we prove that a o n = n o f\inv(Nj)- I^i particular 
this implies that there exists a symbolic dynamics structure on Inv(7V, /). 

The following theorem was proved in [Z3] (see Theorems 5 and 6) for the case n = 2. The 
following is a natural extension to a generic number of sets and the proof is exactly the same. 

Theorem 5.6. The projection tt is onto, and if {xn} G is a periodic sequence, then 
7r~^({2;„}) contains a periodic point. 

We point out that these kind of topological tools have been created to deal with hyperbolic 
periodic points and are not suitable to look for elliptic points. On the other hand the symmetry 
of this system allows the search and proof of (symmetric) periodic points in a much easier way, 
as pointed out in Section 4. Furthermore the main purpose of the paper is to show the existence 
of a symbolic dynamics, which cannot occur in neighborhoods of stable points, therefore we will 
not address this topic any further. 

5.2. Heuristic results. In this section of the paper we provide the heuristic results we obtained, 
while in the following section we prove that such results are rigorous. In the remaining part of 
this paper we fix /i = .3. 
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We need to define 14 t-sets iVo, . . . , iVs, Ni, . . . , N4, Kq, . . . , i^a: the precise definition is given 
in Section 7. Here we only point out that the sets Ni, . . . , N4 are by definition the symmetric 
image with respect to the x-axis of the sets Ni, . . . , N4 (see Definition 5.2), while the remaining 
sets are invariant with respect to the same symmetry. In the following figures the supports of 
the sets A^o — -^5 and their images through the Poincare map are represented, together with a 
portion of the images the supports of the sets Kq and K2. We remark that the images of the 
sets 7V2, Kq and K2 are very thin and appear on the picture as lines. The sets are displayed in 
thick lines, while their images are displayed in thin lines We made two pictures for better clarity. 

Px Px 

0.25r 



0.04 





Lemma 5.7. Let Mi and M2 be t-sets, let Mi and M2 be the sets symmetric with respect to 
symmetry defined above and assume that Mi =^ M2. Then M2 Mi. 

Proof. This follows by the Definitions 5.3, 5.4, 5.2 and the symmetry of the Poincare map (see 
the end of Section 3). □ 

In the following we denote by =^ (resp. ■^=) the covering (resp. the backcovering) with 
respect to the map P and we apply the topological theorems introduced above to the map P. 
The numerical experiments suggest that the following covering relations hold: 



NQ^NQ^Ni^N2^Ns^Ni^N5^ TVg, 



(5.1) 



Kq 
Ki 
K3 



Kq 

K2 
K. 



N3, 
^3, 



If these relations could be verified, then by Lemma 5.7 the following covering relations would 
hold as well: 



(5.2) 



N2^ Ni^ Nq, 

7V3 ^ Kq, 

N3^K2^Ku 
Ns ^ Ks. 
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5.3. Rigorous covering relations. As we pointed out before, it is not convenient to check 
directly the covering relations given in the previous section, since that would require an enormous 
amount of computer time. This problem is due to the fact that, even if we start from a very small 
interval in the computation of the Poincare map, the error and the wrapping effect accumulated 
in the Poincare time is quite large. Instead we prefer to define 12 auxiliary t-sets Mj, i = 0, . . . , 6 
and Lj, « = 1, . . . , 5 on the symmetric Poincare plane, i.e. the plane y = 0, y < 0, and prove 
the covering relations given below (the precise definition of the t-sets is given in Section 7): 

But even to prove these covering relations with the half Poincare maps turns out to be a 
difficult task, indeed some of the sets are quite close to the primaries (particularly the sets L4 
and L5), and this causes the computation of rigorous bounds for the intersection of the trajectory 
with the Poincare plane to be critical. To overcome this difficulty we have to use the alternative, 
but equivalent, definition of covering through the inverse of the Poincare map or backcovering, 
see Definition 5.4. The covering relations we actually prove by computer assistance are as follows: 

TVo ^ Ml ^ TVi ^ M2 ^ 7V2 ^ M3 ^ 7V3 ^ M4 ^ 7V4 ^ M5 ^ TVs 
(5.3) No^Mo, Me^N5, Ko ^ L3 ^ N3 

Note that the coverings L3 4^ Kq, <^ L4, Mq <^ A^o and A^s <^ Mg follow from the 

coverings Kq L3, L4 K^, Nq Mq and Me N5 respectively, by the symmetry of 
the t-sets and Lemma 5.7. As a result, we have to check 21 covering relations. 

In order to have covering relations as described in Definition 5.3, we would need to check 
the image through the Poincare map of the whole support of the t-sets. On the other hand, 
by Lemma 5.1 it is enough to perform rigorous computations on the boundaries of the t-sets, 
provided that we can prove that the map is defined on the whole supports. This is essential for 
the rigorous proof to be made in a reasonable time, since to check by computer assistance the 
definition of the Poincare map on the supports of the t-sets would be very time consuming. In 
order to prove that the half Poincare maps are defined on the supports of all t-sets we argue as 
follows. We only give the proof for the half Poincare map Hi, the proof for H2 being equivalent. 
First we have to prove that the trajectory does not collide with one primary. Then we observe 
that the projections of the trajectories we are interested in on the {y,Py) plane appear to be 
rotating around either a primary or the Lagrangian point Li. If we can compute the angular 
velocity in the (y-Py) plane and prove that it is bounded away from for long enough time, then 
it follows that the trajectory has to cross the Poincare section eventually, hence the map Hi is 
well defined. In [AZ] it was possible to prove by analytical computations that the projection 
on the {x,Px) plane of all trajectories have positive angular velocity with respect to the origin. 
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For the restricted 3-body problem such a general statement does not hold, therefore we have to 
proceed with a different method. First we compute the trajectory of the boundary of the t-set 
with rigorous bounds on the error. The trajectory of the whole boundary describes a "tube" 
in W^. More precisely, the intersection of the trajectory of the boundary with any liyperplane 
(^Py + /^y = bounds a region on the energy surface; the union of all these intersections is what 
we call the tube. Since no trajectories can intersect, all trajectories starting in the interior of 
the support of the t-set have to remain inside the tube. If we can prove that the angular velocity 
in the interior of the tube is always bounded away from 0, say its absolute value is larger than 
e > 0, it follows that the half Poincare map is well defined on the whole interior of the support 
of the t-set; indeed any trajectory starting in the interior of the t-set must intersect the Poincare 
plane in a time T < n/e. 

For a given point {x,px) G dN let 

Mx,Px) = <^(a;,p^;[0,Ti]), 

where is the map defined in Section 3 and Ti = Ti{x,px) is the half Poincare time also 
defined in Section 3. ^{x,Px) is the trajectory of the point {x,Px, 0, ^/2n{x,0) +h-pl) under 
the flow induced by the equations from time to the time Ti when it reaches the Poincare 
plane. If we can prove that Hi{dN) exists, since the flows of two different points cannot 
intersect, unless they coincide, then Hi{dN) bounds a region M in the Poincare plane. We 
want to prove that Hi is defined on the whole set \N\ and Hi{\N\) = M. Let S = \N\ U 
(J <^{x,pg;) U M. The hypersurface E C M"^ divides ^ into two connected regions; call S 
(x,p,)eaiv 

the bounded region. It is clear that for all {x^p^;) G |7V| \ dN either there exists Hi{x,px), or 
(y3 (^2;,Pj;,0, y2?](x7oyT7jr^;T^ G S for all T > 0. Let I3{x,px,y,py) = {2px + ny{x,y))y - 

Py. The next Lemma shows that, if S does not intersects the plane {y = 0, py = 0) and 
\P{x,Px,y,Py)\ > 5 > for all {x,Px,y,Py) € S, then the second case cannot happen. 

Lemma 5.8. Fix h and {x,px) such that x / {it!i,— i?2} o,nd 2Q{x,0) + h — p'^ > 0; let 
Py = ^/2i}(x, 0) + h — p|. Let (x(T),px{t),y{T),py(T)) = ip{x,px,0,Py]T) he the flow induced 
by (1.1); let (3{T) = {2px{T) + ^y{x{T),y{T)))y{T) - pl{T). If there exists ^ > such that 
y'^{T) +pI{T) > and |/3(T)| > 5 for all T G [0,Tr/6] (resp. T G [-7r/6,0]), then the half 
Poincare map Hi (resp. H2) in {x,Px) is well defined. 

Proof. Since y'^{T) + Py{T) > we can use polar coordinates. Let a{T) be the angle. Taking 
the derivative we have 

PyV - pI (2px + ^y)y - pI 

Oi = — = —. 

Since y^ +Py is bounded and {2px + ^y)y —Py is bounded away from 0, then for some T we have 
\a{T) — a{0)\ > TT and the map is defined. □ 

Lemma 5.9. For all covering relations listed in (5.3) the half Poincare maps are defined on all 
the supports of all t-sets. 

Proof. Consider a covering relation with the Hi map, the H2 map and the back-covering relations 
being equivalent. The proof is by computer assistance and it is performed with the following 
procedure. Choose a t-set and let {x,Px) the middle point of its support. Let (p{x,Px,t) be 
the approximate flow, computed by some suitable algorithm (we used Mathematica 4.0 for this 
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purpose). Let Ti be the approximate Poincare time, i.e. (p3{x,Px,t) > for all < i < Ti and 
'f3{x,Px,Ti) = 0. For an integer iV > let = ifi/N, i = 0,N, 

Ei = {{x,px,y,py) G M'^ : d{{x,px,y,Py),(p{x,px,ti)) < 1} , 

where d is a distance in defined by 

j/t 1 1 1 1\ r 2 2 2 2\\ /- I 1 2| i 1 2| i 1 2| i 1 2|\ 

d{{x ,Px,y ,Py),{x ,Px,y ,^2/)) = ™ax(ai|a; -x \,a2\Px - Px\,a3\y -y \,a4\Py-Py\) 
and {flj} are positive constants. Let 

g= U 

i=0,N 

S is a very rough approximation of the flow of the support of the t-set. If all the constants 
entering in the definition of S are chosen appropriately, one can expect that the true trajectory 
of the support of the t-set is entirely contained in S. With computer assistance we checked 
rigorously that the true trajectory of the boundary of the t-sets never leaves S. Using interval 
arithmetics algorithms we compute 

Si = max _ {2px + ^y{x, y))y - p^. 

By computer assistance we prove that maxj 6i is strictly negative for all t-sets. The proof is 
complete by Lemma 5.8. □ 

The previous lemmas yield the following: 
Theorem 5.10. All covering relations listed in (5.3) hold. 

Proof. We check with computer assistance that conditions a' and b in Lemma 5.1 or conditions 
c and d in Definition 5.4 are verified for all the covering relations listed in (5.3). The proof is 
completed by Lemma 5.9. □ 

Consider the collection of the t-sets Nq, Ni, N2, -/V3, JV4, 7V5, Ni, N2, N^, Kq, Ki, K2, 
in this order and let 



T = 



■ 1 


1 



































■ 








1 












































1 












































1 












































1 









































1 











1 














1 



























































1 












































1 








1 





1 


1 


























1 


























1 




















1 















































1 














1 























1 








. 








1 





























1 . 



Theorem 5.10 and the definition of generic covering yield the following 



18 



G. ARIOLI 



Corollary 5.11. All relations in (5.1) and (5.1) hold with generic covering instead of covering 
or backcovering and T is the transition matrix associated with the map P. 

5.4. Symbolic dynamics. Let T^a be the set of all admissible sequences of fourteen symbols 
{No, Ni, N2, Ns, N^, N5, Ni, N2, Ns, N^, Kq, Ki, K2, Ks) with respect to T (recall Definition 
5.8). Corollary 5.11 yields the following result on symbolic dynamics: 

Theorem 5.12. For any biinfinite sequence {a;„} G T,a there exists an orbit of (1.1) which 
crosses the Poincare plane in the t-sets Nq, . . . , Nt^, Ni, . . . and Kq, . . . , Ky, in the order 
prescribed by the sequence. Furthermore, if {xn\ is periodic of period m, such orbit is periodic 
as well and its trajectory on the Poincare section also has period m. 

Proof. The supports of the t-sets Nq, . . . , iVi, . . . , 7V4 and Kq, . . . , K3 are disjoint; let TV 
be the union of all such supports and let tt : Inv(A'^, P) — ^ Eyi be the projection defined as in 
Definition 5.9. The result follows by Theorems 5.6 and Corollary 5.11. □ 

Corollary 5.13. Equation (1.1) admits infinitely many periodic solutions. 

This result is important since it gives the complete picture of the chaotic behavior of the 
system and it yields the result on topological entropy we present in the next subsection. On the 
other hand, a better visual image of the result can be obtained by the following observations. 
Recall that the sets A^o, -^5, Kq and Ks P— cover themselves and Ki P^-covers itself. Hence 
there exists (at least) a periodic orbit for each of these sets which crosses the hyperplane y = 
with positive y velocity in the set (or alternatively in the sets Ki and K2). The orbit crossing 
Nq is the orbit symmetric to U2 with respect to the y axis, while the orbits crossing TVs, Ki 
and are U2, L, Di and Ui respectively. We point out that without considering the derivative 
of the Poincare map we cannot exclude that other periodic orbits cross those t-sets. It would 
be indeed rather easy to rigorously compute the derivative of the Poincare map, but this is not 
relevant for the discussion that follows and such computation will be treated in another paper. 

From the transition matrix it also follows that A^o ^5 ^^=^ -^0 and K-i ■^=^ Nt^ ■^=^ Ki, 
i = 0,2,3, therefore if we consider the t-sets Nq, A^5, Kq, K2 and iCa the transition matrix 
with respect to the map P^^ is a 5 x 5 matrix with all entries equal to 1, that is each of these 
t-sets P^°— covers itself and all the other t-sets. Using again Theorem 5.6 we infer that for all 
biinfinite sequences of five symbols there exists a solution of equation (1.1) which crosses the 
t-sets in the same order. In other words, the system admits solutions which come close to any of 
these orbits in any prescribed order. We point out that this result concern the thirtieth interate 
of the Poincare map, therefore only one crossing of the Poincare section every 30 should be 
considered. 

Finally, we point out that all these results concern the existence of orbits, not uniqueness. 

5.5. Topological entropy estimates. The results presented in the previous subsection yield 
a lower estimate of the topological entropy of the map P by the following lemma: 

Lemma 5.14. Let f : X ^ X be a continuous map. Let S C X be an invariant set, let A be 
an n X n matrix such that there exists a surjective map n : S ^ satisfying a o n = n o f . 
Then the topological entropy of f is larger than ln(max{|A2|, Aj is an eigenvalue 0/^4}). 

Proof. The proof is an easy consequence of Theorem 7.13 in [WA]. □ 

Corollary 5.15. The topological entropy of P is larger than 1.62746. 
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Proof. By Lemma 5.14 a lower bound for ht is given by the maximal norm eigenvalue of the 
transition matrix T. The characteristic polynomial is p{x) = 1 — 2x + 2x^+x^ — "ix^ — 'ix^ + 
7rE'^-4a;*^ + 4a;^-5a;i° + 5a;i2_43.i3_^^i4^ since p(l. 62747) > and p(l. 62746) < 0, then 

ht{P) > 1.62746. □ 

6. Further developments 

Although the 3-body problem is very old, there are still many open problems. Computer 
assisted techniques can give rigorous proofs to some conjectures. This paper is focused on 
results which require only computations. We think that a much richer symbolic dynamics 
can be proved with these methods at different energy levels: this paper should only be considered 
as an example, although in the author opinion Theorem 5.12 and Corollary 5.15 have a value 
in themselves. The rigorous computation of the derivatives of the Poincare map can provide 
information on the linear stability of periodic orbits and yield existence results for homoclinic 
and heteroclinic orbits. We will treat these topics in a forthcoming paper, together with the 
dependence of the results presented in this paper to the ratio of the masses of the primaries. 
A further development under study concerns the study of bifurcations with the energy as a 
parameter value. 

7. Computational details 

7.1. Description of the t-sets. The procedure used to build all t-sets is described in [AZ] for 
the Henon-Heiles Hamiltonian and we refer to that paper for the details. Here we only mention 
that the t-sets Ni and Ki are either centered on the fixed or period 2 points of the Poincare map, 
or placed along the invariant manifolds of such points with the sides approximately parallel to 
the manifolds. The sets Mj and instead are built "by hand" trying to interpolate the image of 
the other sets. We recall that the Poincare section is not connected, since the lines x = ±2^^/^ 
corresponding to the position of the primaries must be excluded. It follows that the invariant 
manifolds of points belonging to different connected components cannot cross. This fact does 
not influence the topological method employed here, since a t-set and its image may be on 
different connected components, bypassing the fact that the Poincare section passes through the 
primaries. The actual t-sets used in the proofs are as follows: 

Definition 7.1. The triple sets are defined by giving the coordinates of the center {x,y), the 
length of the sides {lx,ly) cmd the angular coefficients of the sides (a, /3) as follows: 

Nq: {x,y) = (-0.7942,0.), {l^Jy) = (.007, .007), {a,/3) = (.9325, -.9325). 

Ni: {x,y) = (-.7849, .0129), {la,,ly) = (.002, .005), (a,^) = (-2.136,2.136). 

N2: ix,y) = (-.7589, .05086), {l^.ly) = (.002, .005), (a,/3) = (-1.812,1.981). 

N3: {x,y) = (.08,. 125), ikjy) = (.01, .04), («,/3) = (1.124,2.009). 

N4: {x,ij) = (.131,. 028), (l^^Jy) = (.0076, .012), {a,f3) = (-2.09,2.11). 

N5: {x,y) = (.1471,0.), (l^Jy) = (.015, .015), {a, (3) = (-2.07,2.07). 

Kq: (x,y) = (.04873,0.), (kJy) = (.0002, .0002), (a,/3) = (1.818,-1.818). 

Ki: (x,y) = (.04775,0.), (l^Jy) = (.0009, .0009), (a,/3) = (1.131,-1.131). 

K2: (x,y) = (-.7365,0.), (l^Jy) = (.0012, .0012), («,/3) = (-2,2). 

K3: {x,y) = (.06712,0.), (l^Jy) = (.01,. 01), («,/3) = (1.835,-1.835). 

Mo: {x,y) = (-.1471,0.), (l^^Jy) = (.0065, .0065), (a,/3) = (1,-1). 

Ml.- {x,y) = (-.1424, .007289), (l^Jy) = (.0026, .0065), {a, (3) = (.9,-1.05). 

M2: {x,y) = (-.1311,. 0279), {l^,ly) = (.0025, .005), (a,^) = (4.2,2). 
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M3; {x,y) = (-.0792,. 1267), {l^.ly) = (.002, .005), («,/3) = (4.2,1.8). 

M4.- (a;,y) = (.7592, .05021), {l^,ly) = (.006, .018), (a,/3) = (1.2,2.2). 

Ms.- {x,y) = (.7854, .01406), {l^,ly) = (.01,. 008), (a,/3) = (4.2,2.2). 

Me.- {x,y) = (0.7942,0.), {l^Jy) = (.013, .013), («,/3) = (1,-1). 

Li: {x,y) = (-.0742, -.1015), {l^,ly) = (.0009, .0009), (q;,/3) = (1.131,-1.131). 

L2: {x,y) = (-0.0738, .1025), {l^,ly) = (.002, .002), (a,^) = (-1.85,1.85). 

L3.- {x,y) = (-.04873,0.), (l^Jy) = (.01,. 01), (a,/3) = (1.819,-1.819). 

L4.- ix,y) = (.6563,0.), (l^Jy) = (.005, .005), («,/3) = (-1.49,1.49). 

L5.- {x,y) = (.6545, .022), (l^Jy) = (.01,. 01), (a,/3) = (-1.49,1.49). 

The left (resp. right) edge of each set is the segments whose end points are (x + l^cosa + 
lyCOS /3,y + lx sin a + ly sin /3) and (x + lx cos a — ly cos /3,y + lx sin a — ly sin /3) (resp. (x — lx cosa + 
ly COS f3,y — Ix sino; + ly sin^) and {x — Ix cos a — ly cos f3,y — Ix sin a — ly sin/3) ). The boundaries 
of the left and right sides of each t-set are the lines crossing two opposite vertices of the support, 
see Figure 3. 

Furthermore, we denote Ni, i = 1, ... 4, the sets obtained by reflecting the sets N-i with respect 
to the x-axis, see Definition 5.2. Note that the sets Nq, N5, Ko, Ki, K2, K3, are by definition 
symmetric with respect to the x axis. 

7.2. Computation techniques. We describe here the algorithm used in the computer assisted 
proof of Theorem 5.10. The proof of Theorem 4.1 is equivalent. By Lemma 5.1 the proof consists 
in checking that the images through the Poincare map of the edges of the t-sets lie in some 
assigned regions of the plane as described by the covering relations. We do not know the exact 
images of such edges, since no analytical solution of the equation is available, therefore all we 
can do is to estimate the trajectories and compute the intersections with the Poincare plane 
with rigorous error bounds. In order to compute the image of a side, we partition it in segments 
with small enough length and we check that every such segment is mapped in the correct region. 
To compute the image of a small segment, we enclose it in an interval set (a rectangle) and 
we compute its trajectory with a Taylor-Lohner algorithm using interval arithmetics. More 
precisely, we start with a Taylor method of order 12, i.e. we estimate the trajectory of an 
interval by using the Taylor expansion of order 12 and we estimate the error by the Lagrange 
remainder. If h is the time step, we compute a rough but rigorous enclosure D of the trajectory 
at times [0, /i], that is an interval set D such that the solution of the equation lies in D for all 
times between and /t, and by Lagrange theorem we estimate the error we make neglecting 
the remaining terms of the Taylor expansion by computing x^^'^\D)'-j^ where x'^^'^\D) (which 
is an interval enclosing all possible values assumed by the 13th derivative of the trajectory, 
therefore enclosing the Lagrange remainder) is computed using a recursive algorithm for the 
time derivatives of the solutions. 

The interval arithmetics algorithms address the problem of computing the trajectory of a 
segment and of keeping track of the errors in an elegant and rigorous way, but they introduce 
another problem. Indeed, even in the simplest dynamical system, the procedure described above 
leads to a very rough estimate of trajectories, due to the wrapping effect which makes the bounds 
on the error grow exponentially. The problem has been strongly contained by introducing the 
half Poincare maps and the backcovering relation, see the discussion in Section 5, but these 
techniques do not suffice. We obtain another significant reduction of the wrapping effect by 
using the Lohner algorithm. We refer to see Section 6 in [AZ] and references cited there for 
a discussion of interval arithmetics and wrapping effect and for a description of the Lohner 
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algorithm employed. See also [MZ] for a discussion on interval arithmetics and [L] for the Lohner 
algorithm. We point out that using the Lohner algorithm to compute directly the Poincare map, 
instead of the half Poincare maps or their inverse, the computation time to perform all the proofs 
is unrealistic on current desk computers, therefore the definitions of the half Poincare maps and 
of backcovering are essential to perform the proofs. 

The round-off errors are taken care directly by suitable C++ libraries and by Mathematica. 
Such errors may vary by changing computer and/or operationg system, but since they are usually 
very small when compared to the wrapping effect and since all proofs go through with a relatively 
large margin, we expect that the proof can be easily reproduced on any recent computer. 

The typical time step used in the computation of the images of the t-sets is dt = 10~^, but 
in few cases we had to use some lower value, down to dt = 5.10"^. Each side of the t-sets has 
been usually divided in 100 to 2000 segments, depending on the apparent value of the Lipschitz 
constant of the map in the area considered. In a few difficult cases we had to divide each side of 
a t-set in 5000 segments. Most of the computations used to obtain the bounds for the location 
of the orbits in Table 1 used dt = 2.10"^. 

To perform the proofs the author implemented a version of the whole algorithms in a combi- 
nation of Mathematica and C++ under the Linux O.S. More precisely, Mathematica has been 
only used to handle all the data and to perform a few algorithms which are less demanding 
for the CPU, but more complicated to implement. Furthermore Mathematica has been used 
to make all numerical experiments used to choose the t-sets and to draw the pictures. On 
the other hand C++ has been used for the heavy interval arithmetic computations, where it 
offered a much better speed. The connection between the two languages is obtained by Math- 
Link. We wish to point out that the full proof took almost a month of CPU time on a machine 
equipped with a IGHz Pentium III processor. In fact a large amount of the time is used for the 
computations involving the set -L5 whose trajectory is very close to one primary. A full C++ 
algorithm would reduce the time at the price of a much more complicated and less user-friendly 
programming. We think that this sharing of tasks is almost optimal, as far as computational 
speed and simplicity of programming and data handling is concerned. The reader who desires 
to reproduce the computer assisted proofs in this paper without writing the program can use 
the Mathematica notebook which is attached to this preprint. The notebook is provided with 
comments and instructions. By using the commands in the notebook it is possible to make an 
independent computation of the images of the sides of the t-sets and an independent check of 
Theorem 4.1. 

Acknowledgement. The author is very grateful to P. Zgliczynski for many discussions and 
for the interval arithmetic C++ libraries. 
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